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LONG-TERM  GOAL 

Our  long  term  goals  are  to  improve  the  modeling  and  prediction  capabilities  of  numerical  ocean  models 
and  data  assimilation  systems  in  support  of  operational  coastal  and  ocean  forecasting. 

OBJECTIVES 

We  wish  to  systematically  evaluate  the  skill  of  basin-scale  and  coastal  ocean  models,  in  producing  the 
observed  climatologies  and  variabilities  of  the  North  Atlantic  Ocean.  We  also  wish  to  test  data 
assimilation  methodologies,  in  particular  those  involving  satellite-derived  surface  data,  in  order  to 
improve  the  skill  of  nowcasting  and  forecasting  meso-scale  features  in  the  open  ocean,  in  the  Gulf 
Stream  region  and  the  coastal  zone. 

APPROACH 

For  modeling  and  data  assimilation  studies  we  rely  on  the  sigma  coordinate  Princeton  Ocean  Model, 
POM  (Blumberg  and  Mellor,  1987;  Mellor,  1996)  using  various  configurations  ranging  from  high 
resolution  coastal  ocean  grids  to  a  low  resolution  climate  grid  of  the  entire  Atlantic  Ocean.  The 
approach  of  using  different  model  domains  allows  us  to  inmprove  open  boundary  conditions  in  regional 
models  by  using  information  obtained  from  larger  domains.  Curvilinear  orthogonal  grids  are  used  to 
provide  higher  resolutions  near  the  coast  and  near  eddy-rich  regions.  Sensitivity  studies  evaluate  the 
sensitivity  of  the  model  climatology  and  variability  to  various  parameters  and  configurations.  The  data 
assimilation  approach  is  based  on  previously  developed  optimal  interpolation  methodologies  of  using 
satellite-derived  surface  data  to  update  subsurface  fields  (Mellor  and  Ezer,  1991;  Ezer  and  Mellor, 
1994).  The  emphasis  during  the  reported  period  was  on  two  new  areas  of  research:  1.  To  test  the 
assimilation  approach,  previously  tested  only  with  regional  Gulf  Stream  models,  with  basin-scale  models 
and  2.  To  test  the  approach  within  a  real-time  operational  coastal  forecast  system.  An  important  part  of 
the  research  was  done  through  our  participation  in  the  Coastal  Ocean  Forecast  System,  COFS,  project 
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(Aikman  et  al.,  1996)  and  in  the  Data  Assimilation  and  Model  Evaluation  Experiments  -  the  North 
Atlantic  Basin  (DAMEE-NAB). 

WORK  COMPLETED 

A  basin-scale  North  Atlantic  model,  extending  from  5°N  to  50°N,  has  been  set  up  and  tested,  as  part  of 
an  intercomparison  study  (DAMEE)  that  involved  two  sigma  coordinate  models,  two  z-level  models 
and  two  layer  models.  We  evaluated  the  sensitivity  of  POM  to  various  parameters  and  configurations 
and  compared  the  results  with  observations  and  with  other  models.  The  lateral  boundary  condition  plays 
an  important  role  in  the  model  climatology  and  variability;  thus  we  compared  closed  and  open  boundary 
conditions  and  found  that  better  Gulf  Stream  separation  and  more  realistic  recirculation  gyres  are 
obtained  with  open  boundary  conditions  using  transports  imposed  from  a  larger  model  domain 
compared  with  closed  boundary  conditions  and  a  buffer  zone,  a  result  consistent  with  previous  Gulf 
Stream  separation  studies  using  a  smaller  domain  (Ezer  and  Mellor,  1992).  Two  different  grid 
resolutions  were  tested,  both  having  a  curvilinear  orthogonal  grids  with  higher  resolution  in  the  Gulf 
Stream  and  the  Gulf  of  Mexico  regions  and  lower  resolution  grids  in  the  eastern  Atlantic  region.  Since 
POM  uses  the  grid-dependent  Smagorinsky  diffusivity  formulations,  each  grid  has  been  tested  with 
different  values  of  diffusion  and  viscosity,  thus  with  different  values  of  the  turbulent  Prandtl  number.  Of 
particular  importance  to  sigma  coordinate  models  is  the  evaluation  we  performed  on  the  effect  of  the 
along-sigma  diffusion  formulation  on  the  model  climatology  (Ezer,  1998). 

Our  data  assimilation  scheme  has  been  tested  with  a  basin-scale  POM  domain  for  the  first  time, 
assimilating  three  years  of  Topex/Poseidon  altimeter  data,  from  1993  to  1995,  into  the  North  Atlantic 
ocean  model  and  evaluating  the  assimilated  model  against  unassimilated  model  and  against  observed 
fields  (Ezer  and  Mellor,  1998).  A  procedure  to  assimilate  frontal  position  data,  such  as  the  Gulf  Stream 
north  wall  analysis  derived  from  AVHRR  images,  has  been  completed  and  soon  will  be  integrated  within 
the  operational  forecast  system  (COFS)  running  at  NCEP,  using  the  U.S  Navy’s  daily  frontal  analysis 
data. 

The  problem  of  the  pressure  gradient  errors  in  sigma  coordinate  models  (Mellor  et  al.,  1994)  has  been 
revisited  and  analyzed  in  more  detail.  The  study  shows  that  previous  descriptions  of  these  errors  are 
incomplete  or  are  in  error  (Mellor  et  al.,  1998a).  Development  of  a  generalized  vertical  coordinate 
configuration  of  POM  has  been  completed  (Mellor  et  al.,  1998b);  it  will  enable  the  combined  use  and 
the  intercomparison  of  z-level,  sigma-level  and  possibly  isopycnal  models  using  otherwise  identical 
numerics. 

RESULTS 

Sensitivity  studies  with  a  basin-scale  POM  configuration  show  that  the  effect  of  the  Smagorinsky 
horizontal  diffusivity  and  viscosity  on  the  model  variability  becomes  less  significant  as  resolution 
improves.  A  change  in  grid  size  of  50%  has  a  larger  effect  on  the  variability  than  an  order  of  magnitude 
change  in  the  diffusivity.  The  along-sigma  diffusion  together  with  removal  of  climatological  fields 
results  in  a  slight  reduction  of  the  model  long-term  climate  drift,  but  does  not  affect  the  short-term 
seasonal  and  mesoscale  variabilities;  in  fact,  the  results  of  a  standard  Smagorinsky  parameterization  are 
very  similar  to  those  of  no  diffusion  at  all  (the  sigma  coordinate  model  was  found  to  be  numerically 
stable  even  with  zero  horizontal  diffusion).  The  model  variability  compares  well  with  the  observed 
variability  only  in  areas  where  model  grid  is  small  enough  (around  10  km)  such  as  in  the  Gulf  of  Mexico 
and  in  portions  of  the  Gulf  Stream. 
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Fig.  1.  Surface  elevation  over  the  Gulf  of  Mexico  and  the  Caribbean  Sea  on  10  August,  1994  from 
basin-scale  North  Atlantic  model,  (a)  Model  prediction  without  assimilation,  (b)  model  hindcast 
with  assimilation,  and  (c)  Topex/Poseidon  altimeter  data  interpolated  into  the  model  grid  without 
involvement  of  model  dynamics.  Note  that  the  assimilated  model  produces  the  observed  meso- 
scale  features  while  maintaining  a  more  realistic  Florida  Current  than  that  of  the  altimeter  analysis. 


Data  assimilation  studies  show  considerable  improvement  in  the  model  variability  compared  to  an 
unassimilated  model  and  some  nowcast  skill  in  producing  the  observed  mesoscale  features,  but  only  in 
regions  with  sufficient  model  resolution.  The  Topex/Poseidon  altimeter  data  were  used  together  with 
surface- subsurface  correlations  to  continuously  update  the  temperature  fields  in  the  entire  water  column 
and  then  the  model  surface  elevation  was  compared  with  objective  analysis  of  the  altimeter  data  without 
a  dynamic  model  in  order  to  test  the  vertical  projection  approach  and  the  contribution  of  the  dynamical 
model.  An  example  of  such  a  comparison  is  shown  in  Fig.  1,  for  a  small  portion  of  the  North  Atlantic 
model  domain;  this  example  shows  the  importance  of  the  model  dynamics  in  improving  the  nowcast 
based  on  data  alone. 

IMPACT/APPLICATION 

The  improvements  in  the  numerical  model  itself  and  in  the  data  assimilation  techniques  will  continue  to 
advantagely  affect  the  POM  user  community  in  particular,  and  prediction  systems  in  general.  For 
example,  the  now  accepted  fact  that  one  must  project  surface  information  into  the  entire  water  column 
followed  our  1991  paper  on  the  subject.  The  assimilation  scheme  that  is  soon  to  be  implemented  in  the 
operational  COFS  will  provide  an  important  test  bed  for  other  Navy’s  forecast  systems  that  use  POM. 

TRANSITIONS 

The  transition  of  model  and  data  assimilation  algorithms  as  well  as  knowledge  from  the  research  and 
development  stage  to  the  operational  environment  has  been  going  on  during  the  past  few  years  and  will 
continue.  Continuous  transition  of  model  developments  from  the  Princeton  group  to  the  ocean  modeling 
community  is  done  through  our  internet-based  POM  users  group,  which  includes  more  than  360 
registered  users  from  30  countries. 

RELATED  PROJECTS 

We  have  fruitful  collaborations  with  many  scientists  through  our  involvement  in  the  National  Ocean 
Partnership  Program  (NOPP),  the  Data  Assimilation  and  Model  Evaluation  Experiments  (DAMEE),  the 
Coastal  Ocean  Forecast  System  (COFS)  for  the  U.S.  east  coast  and  the  POM -based  forecast  system  for 
the  U.S.  west  coast.. 
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